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Pollination is an early and critical step in plant reproduction, 
leading to successful fertilization. It consists of many sequen- 
tial processes, including adhesion of pollen grains onto the 
surface of stigmatic papilla cells, foot formation to 
strengthen pollen-stigma interaction, pollen hydration 
and germination, and pollen tube elongation and penetra- 
tion. We have focused on an examination of the expressed 
genes in papilla cells, to increase understanding of the mo- 
lecular systems of pollination. From three representative 
species of Brassicaceae (Arabidopsis thaliana, A. halleri and 
Brassica rapa), stigmatic papilla cells were isolated precisely 
by laser microdissection, and cell type-specific gene expres- 
sion in papilla cells was determined by RNA sequencing. As a 
result, 17,240, 19,260 and 21,026 unigenes were defined in 
papilla cells of A. thalicma, A. halleri and B. rapa, respectively, 
and, among these, 12,311 genes were common to all three 
species. Among thel 7,240 genes predicted in A. thaliana, 
one-third were papilla specific while approximately half 
of the genes were detected in all tissues examined. 
Bioinformatics analysis revealed that genes related to a 
wide range of reproduction and development functions 
are expressed in papilla cells, particularly metabolism, tran- 
scription and membrane-mediated information exchange. 
These results reflect the conserved features of general cellu- 
lar function and also the specific reproductive role of papilla 
cells, highlighting a complex cellular system regulated by a 
diverse range of molecules in these cells. This study provides 



fundamental biological knowledge to dissect the molecular 
mechanisms of pollination in papilla cells and will shed light 
on our understanding of plant reproduction mechanisms. 

Keywords: Arabidopsis halleri • Arabidopsis thaliana • Brassica 
rapa • Laser microdissection • Papilla cell • RNA sequencing. 

Abbreviations: bHLH, basic helix-loop-helix; BR, brassinos- 
teroid; GO, gene ontology; LM, laser microdissection; LRR- 
RLK, leucine-rich repeat receptor-like protein kinase; PME, 
pectin metylesterase; qRT-PCR, quantitative real-time PCR; 
RIN, RNA integrity number; RNA-seq, RNA sequencing; SI, 
self-incompatibility; TF, transcription factor. 



Introduction 




In flowering plants, several steps occur from anthesis to fertil- 
ization during sexual reproduction. The fertilization process 
begins when the pollen grains land on the stigma surface, 
that is the papilla cells. For pollen to adhere to the stigma, 
proteins and lipids are provided from the stigma and pollen; 
these mix to form a structure called the foot (Gaude and 
Dumas 1984, Dickinson 1995). In particular, proteinaceous 
components and lipids of the pollen coat are essential for 
both foot formation and subsequent pollen hydration, 
termed 'coat conversion' (Elleman and Dickinson 1986). After 
foot formation between pollen and papilla cells, water and 



Plant Cell Physiol. 54(11): 1894-1906 (2013) doi:10.1093/pcp/pct133, available online at www.pcp.oxfordjournals.org 

© The Author 2013. Published by Oxford University Press on behalf of Japanese Society of Plant Physiologists. 

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http:// 

creativecommons.org/licenses/by-nc/3.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, 

provided the original work is properly cited. For commercial re-use, please contactjournals.permissions@oup.com 



1894 Plant Cell Physiol. 54(11): 1894-1906 (2013) doi:10.1093/pcp/pct133 © The Author 2013. 



Papilla-specific transcriptome in Brassicaceae 

PLANT & CELL PHYSIOLOGY 



resources for germination and pollen tube elongation are pro- 
vided from the stigma to pollen, through a capillary system 
(Zuberi and Dickinson 1985, Zinkl et al. 1999). Once the 
pollen grain has germinated, the pollen tube directionally pene- 
trates into a papilla cell, elongates through the stigma and style, 
and enters into the ovary (Cheung et al. 1995, Franklin-Tong 
1999, Lord 2000, Shimizu and Okada 2000). The pollen tube 
eventually arrives at the embryo sac in the ovule, and finally, in 
the ovary, two sperm cells are released from the pollen tube tip 
and fertilize an egg cell and central cell, respectively, resulting in 
seed development (Dickinson and Elleman 1994, Stephenson 
et al. 1997). 

In these plant reproduction processes, pollination is an 
early and important step for successful fertilization. Papilla 
cells have various roles in this step, such as discrimination of 
pollen grains for inter- and intraspecies barriers, pollen hydra- 
tion, supply of stigmatic resources for pollen germination 
and pollen tube elongation, and pollen tube guidance (de 
Nettancourt 2001, Goto et al. 2011). However, the molecular 
basis of these mechanisms of pollination is still largely unclear 
(Hiroi et al. 2013). 

Flowering plants have evolved exquisite systems for the con- 
trol of mating interactions before fertilization, and the most 
extensively studied of these is self-incompatibility (SI) (Suwabe 
et al. 2010, Watanabe et al. 2012). SI is a genetic system for 
intraspecific pollen selectivity and is recognized as the most 
sophisticated system to prevent selfing (Bateman 1995). It is 
controlled by a single, highly polymorphic S locus and is mani- 
fested during the interaction of pollen grain with the papilla 
cells of the stigma epidermis (Watanabe et al. 1994, Suzuki et al. 
1999). In Brassicaceae, S-locus receptor kinase (SRK) encodes a 
Ser/Thr receptor kinase, which localizes to the stigmatic plasma 
membrane and interacts with the pollen coat cysteine-rich 
protein (SP11/SCR) to induce the SI reaction (Schopfer et al. 
1999, Takasaki et al. 2000, Takayama et al. 2000, Takayama et al. 
2001). Downstream factors of the SI system, M-locus receptor 
kinase (MLPK), arm repeat containing 1 (ARC1) and Exo70A1, a 
subunit of the exocyst complex, have been identified, which 
mediate pollen-stigma interactions (Gu et al. 1998, Stone et al. 
1999, Stone et al. 2003, Murase et al. 2004, Kakita et al. 2007a, 
Kakita et al. 2007b, Samuel et al. 2009). However, while some of 
these factors have been identified in the SI signaling cascade, 
the overall systems of SI and pollen-stigma interaction at the 
molecular level are still largely unclear, even in the much stu- 
died Brassicaceae. 

For establishment of a knowledge base to aid understanding 
of the molecular mechanisms of pollination in the stigmatic 
papilla cell, including SI, transcriptome analysis is undoubtedly 
informative (Matsuo et al. 2011, Bai et al. 2013). So far, genome- 
wide transcriptome analysis of the stigma has been conducted 
in Arabidopsis, rice and maize, to study pollen-stigma commu- 
nication and pollen tube guidance, by microarray technology 
(Swanson et al. 2004, Li et al. 2007). However, the microarray 
system has several limitations: high background levels due to 
cross-hybridization; lack of sensitivity at low expression levels; 



and reliance upon existing knowledge of genome sequences 
(Huang et al. 2013). In contrast, RNA sequencing (RNA-seq) 
has the potential to overcome these limitations (Wang et al. 
2009) and offers a variety of new possibilities in biology, such as 
whole-genome sequencing, global identification of expressed 
genes, gene expression profiling and identification of novel 
genes, alternative splicing and sequence variation (Marioni 
et al. 2008, Garg et al. 2011). 

Another problem in studying gene expression is that be- 
cause it is difficult to isolate single cell types specifically from 
plant organs, many microarray studies have been conducted 
using RNAs extracted from a mixture of tissues and/or cell 
types (Autran et al. 2011). This limitation makes it difficult to 
analyze expressed genes at high precision in the target cell 
types. To overcome this, laser microdissection (LM) is a valuable 
tool to isolate specific cell types from sectioned specimens of 
heterogeneous tissues (Girke et al. 2000, Asano et al. 2002, 
Nakazono et al. 2003). LM allows the isolation of individual 
cells directly from the surrounding tissue, based on histological 
identification, by laser capture and laser cutting (Kerk et al. 
2003, Suwabe et al. 2008, Kubo et al. 2013). Given the 
advantages and opportunities offered by RNA-seq and LM, 
this combination has the potential for precise and efficient 
analysis of transcriptional gene expression profiles of individual 
cell types. 

Here, we report the establishment and characterization of 
papilla cell expressed gene data sets, among existing gene data 
in public databases, in Brassicaceae species Arabidopsis thali- 
ana, A. halleri and Brassica rapa, by LM- RNA-seq analysis. 
Arabidopsis thaliana is a model plant species and A. halleri is 
a close relative. Brassica rapa, a vegetable member of the 
Brassicaceae, is a model species of the genus Brassica, and the 
annotated draft genome sequence has recently been com- 
pleted (Wang et al. 2011). By using LM technology, the papilla 
cells of these species were isolated precisely, and papilla cell- 
specific gene expression profiles were obtained by RNA-seq. 
Using papilla cell-specific gene expression data sets from A. 
thaliana, A. halleri and B. rapa, we investigated the molecular 
biology of papilla cells by bioinformatics analysis, in terms of 
pollination mechanisms and papilla cell development. 



Results and Discussion 




Optimal conditions for paraffin-embedded 
sectioning of papilla cells 

It is known that high integrity RNA is critical for gene expres- 
sion analysis, and thus we first evaluated the optimal histolo- 
gical fixation conditions for paraffin embedding papilla cells for 
RNA extraction. The use of four fixatives (100% ethanol, 75% 
ethanol/25% acetate, 60% ethanol/40% acetate and 100% acet- 
one) was compared for fixation of stigmas, including papilla 
cells, of A. thaliana, A. halleri and B. rapa. The quality of total 
RNA of LM papilla cells in each fixative was assessed using an 
Agilent 2100 Bioanalyzer. In A. thaliana, RNA quality expressed 
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as the RNA integrity number (RIN) = 7.0, 6.9, 6.2 and 6.5 in 100% 
ethanol, 75% ethanol/25% acetate, 60% ethanol/40% acetate 
and 100% acetone, respectively. In A. halleri, the RIN = 7.1, 
6.8, 7.1 and 6.3 in 100% ethanol, 75% ethanol/25% acetate, 
60% ethanol/40% acetate and 100% acetone, respectively; and 
in B. rapa, the RIN = 5.5, 3.7, 3.5 and 6.2 in 100% ethanol, 75% 
ethanol/25% acetate, 60% ethanol/40% acetate and 100% acet- 
one, respectively (Table 1). Thus, optimal fixation conditions 
varied depending on the plant species. In addition, morpho- 
logical disintegration of tissue sections, especially in the ovary, 
was occasionally observed in 100% ethanol and 100% acetone 
(data not shown). Taken together, the best fixative conditions 
in terms of RNA quality and morphology was 75% ethanol/25% 
acetate in A. thaliana, 60% ethanol/40% acetate in A. halleri and 
100% acetone in B. rapa. 

RNA-seq and identification of papilla-expressed 
genes 

Using optimal fixation conditions for pistils, embedded speci- 
mens were prepared using 2-6 floral buds from each plant 
species. Serial paraffin sections of pistil specimens were pre- 
pared at 6-8 urn in thickness and papilla cells were isolated 
by LM, using approximately 60 pistil sections from each plant 
sample (Fig. 1). Total RNAs were extracted from the LM-col- 
lected papilla cells, confirmed to be acceptable quality for RNA- 
seq analysis (Supplementary Fig. S1), and 1.5 rounds of mRNA 



Table 1 Optimal fixative conditions for different Brassicaceae 
species 



Plant 


RNA quality (RIN) 






species 


100% ethanol 


Ethanol/acetic acid 


100% acetone 






75%/25% 


60%/40% 




A. thaliana 


7.0 


6.9 


6.2 


6.5 


A. halleri 


7.1 


6.8 


7.1 


6.3 


B. rapa 


5.5 


3.7 


3.5 


6.2 



linear amplification were performed for each sample. 
Subsequently, 50-200 ngul -1 cDNAs were obtained from 
each sample, and cDNA fragment libraries were subjected to 
SOLiD5500xl sequencing. As a result of the RNA-seq analysis, 
65,690,858, 101,509,276 and 74,175,637 raw reads of 75 bp in 
length were obtained in A. thaliana, A. halleri and B. rapa, re- 
spectively. On the alignment of raw reads to the A. thaliana and 
B. rapa reference genomes, a total of 21,310, 23,414 and 26,852 
unigenes were defined, respectively, after removing low-quality 
raw reads (MapQ <10) and adaptor sequences. Of these, 
17,240, 19,260 and 21,026 unigenes were found to be >10 
read tags and could be assigned an annotation based on data- 
base sequence assembly in A. thaliana, A. halleri and B. rapa, 
respectively (Table 2; Supplementary Tables S1-S3). In 
A. thaliana, approximately 50-60% of the genes have previously 
been predicted to be expressed in stigmas (Suwanson et al. 
2005), and the remaining ^50% of the genes have never been 
predicted to be expressed in stigmas, and thus they are likely to 
be novel stigma-expressed genes, including papilla-expressed 
genes. The fact that approximately half of the annotated 
genes from the RNA-seq have not been detected previously 
in microarray and cDNA subtraction analyses highlights the 
advantages, including greater sensitivity, of the combination 
of LM and RNA-seq analysis for the tissue of interest. 

Among the top 100 most highly expressed genes of the 
A. thaliana papilla cell transcriptome (Supplementary Table 
S1), gene ontology (GO) annotation analysis showed that genes 
in the categories membrane-related (54/100), cell wall (25/100) 
and metabolism (24/100) were significantly enriched and some 
of these have already been proved to be involved in plant re- 
production systems. For example, among these genes, the SCA- 
like Arabidopsis UP, AT2G38540, has been shown to have vari- 
ous roles both in plant reproduction and in growth (Chae et al. 
2010). PDIL2-1, AT2G47470, which encodes a cell wall-related 
protein disulfide isomerase, has been verified to be important 
in embryo sac maturation (Wang et al. 2008). This demon- 
strates that these categories include genes that are expected 




Fig. 1 Isolation of a papilla cell by LM in Arabidopsis and Brassica. The paraffin-embedded stigmas were cut into 6-8 jam thick sections using a 
microtome and mounted on PEN membrane frame slides. After removing the paraffin, LM was performed using an Arcturus XT Laser Capture 
Microdissection System. (A-C) Arabidopsis thaliana, (D-F) B. rapa, (A, D) before LM, (B, E) after LM, and (C, F) isolated papilla cells. 
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Table 2 Summary of experimental tags sequenced in reference A. thaliana and B. rapa genomes 



Plant species 



Total reads 



Mapped reads 



MapQ >10 



Total genes 



Unigenes 
Reads >0 



Reads >10 



A. thaliana 

A. halleri b 

B. rapa 



65,690,858 

101,509,276 

74,175,637 



38,955,652 
59.30% 

33,613,230 
33.11% 

62,343,265 
84.05% 



35,266,802 
53.69% 

25,431,815 
25.05% 

54,037,867 
72.85% 



33,602 a 



41,173 c 



21,310 
63.42% d 

23,414 
69.68% d 

26,852 
65.22% d 



17,240 
51.31% d 

19,260 
57.32% d 

21,026 
51.07% d 



a Total gene number from A. thaliana deposited in the database (http://www.arabidopsis.org/). 
b The total gene number from A. halleri was not determined. 

c Total gene number from B.rapa deposited in the database (http://brassicadb.org/brad/). 

d Unigenes/total genesx100 (%). For both A. thaliana and A. halleri, the total gene number used for calculation was 33,602. 



A. thaliana 



stigma 
leaf 




AT1G25330 



AT1G65720 



AT1G78040 



AT2G25520 



AT2G33850 



AT3G27810 




Bra019277 



Bra018477 



Bra010085 



Bra028514 Bra012322 



Bra030385 



Fig. 2 Expression analysis of genes identified in the papilla by LM-RNA-seq. qRT-PCR was performed using the SYBR Green Real-Time PCR 
Master Mix with a Bio-Rad CFX96 Real-Time Detection System. cDNAs were prepared from total RNAs of stigma and leaf. Representation of the 
relative expression levels of the genes was normalized to the expression level of UBC21. Quantitative validation of three replicates was calculated 
using the delta-delta threshold cycle relative quantification method. 



to be particularly important in papilla cell processes, such 
as regulation of papilla cell formation and pollen-papilla 
communication. 

Validation of the LM-RNA-seq transcriptome in 
papilla cells 

To assess the validity and reliability of our LM-RNA-seq data 
from papilla cells, quantitative real-time PCR (qRT-PCR) ana- 
lysis was performed on 12 representative genes belonging to 
different classification groups (Supplementary Table S4). 
Expression levels of the 12 genes were compared between pa- 
pilla cells and the leaf to determine whether these genes were 
expressed in papilla specifically or preferentially (Fig. 2). The 



results from qRT-PCR showed that these 12 genes were 
expressed exclusively in papilla cells, and this result correlated 
with those obtained by LM-RNA-seq. To confirm further 
the expression and localization of the genes in papilla cells, in 
situ hybridizations were performed with six representative 
genes, of stress-related, metal tolerance and unknown functions 
(Supplementary Table S5). Probes from all six genes exhibited 
hybridization signals in papilla cells and the transmitting 
tract of the stigma, confirming that the corresponding 
gene transcripts were detected in papilla cells by in situ hybrid- 
ization (Fig. 3). 

Taken together, these results demonstrate the validity 
and reliability of our LM-RNA-seq data in papilla cells. This 
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Fig. 3 Confirmation of gene expression patterns by in situ hybridiza- 
tion. The pistils for in situ hybridization were collected 1 d before 
anthesis. Both antisense and sense probes were synthesized using a 
T7/SP6 digoxigenin RNA labeling kit. Insets show the signal area in 
papilla cells at high magnification. (A) AT1G09690, (B) AT1G53130, 
(C) AT2G41430, (D) Bra004382, (E) Bra019017 and (F) Bra028171. 



Papilla cell 

17,240 




Stamen Leaf 

12,143 9,868 



Fig. 4 Comparison of transcriptome data sets from papilla cells, 
stamen, leaf and flower in A. thaliana. (A) Venn diagram for the 
number of genes expressed in papilla cells, stamen and leaf. 
Expression data for stamen and leaf are from Schmid et al. (2005). 
(B). Venn diagram for the number of genes expressed in papilla cells 
and flowers. Expression data for flowers are from Schmid et al. (2005). 



also highlights the potential capability of LM-RNA-seq as a 
technique for fine-scale gene expression analysis. 

Comparison of expressed genes in papilla cells 
and other tissues in A. thaliana 

To characterize the Arabidopsis papilla-specific genes, LM- 
RNA-seq data from papilla cells were compared with those 
from microarray analysis of stamen and leaf in A. thaliana 
(Schmid et al. 2005) (Fig. 4A). Among the 17,240 genes 
expressed in papilla cells, approximately half (8,538/17,240) 
were uniformly detected in all tissues examined here, suggesting 
that these are conserved genes for general cellular function 
common to all tissues. In contrast, one-third of the genes 
(5,668/17,240) were papilla specific. Because papilla cells have 
specialized functions, focusing on pollination, these genes are 
expected to be involved in the regulation and maintenance of 
the specific roles of papilla cells. 

It is noteworthy that 6,898 papilla cell-expressed genes from 
the present study were not identified in the previous micro- 
array data from whole flowers (Schmid et al. 2005) (Fig. 4B). 
Furthermore, 40% of the genes (2,749/6,898) were also not con- 
tained on the Arabidopsis ATH1 Genome Array GeneChip, 
which includes >22,500 probe sets representing approximately 
24,000 genes. This reflects the limitation of microarray 
techniques, and at the same time highlights the potential of 
LM-RNA-seq technology to increase greatly the accuracy of cell 
type-specific gene isolation and the scale of transcriptome 
analysis. 

Comparison of papilla-expressed gene data sets 
from the three Brassicaceae species 

To reveal the similarities and differences in papilla-expressed 
genes in Brassicaceae, the data sets from A. thaliana, A. halleri 
and B. rapa were compared (Fig. 5). Because the genome of 
Brassica has been triplicated during the course of evolution, 



A. thaliana 

17,240 , 



1,081 



1,607, 



,241 



A. halleri 

19,260 





12,31 r\ 




1 2,376 1 


966 


1 581 J 



B. rapa 

14,099 



Fig. 5 Conservation of papilla-expressed genes in Brassicaceae. Venn 
diagram for the number of genes in papilla cells from A. thaliana, 
A. halleri and B. rapa. 



every gene has essentially three copies, i.e. one orthologous 
and two paralogous genes, in the genome. Thus, in the case 
of B. rapa, 14,099 unigenes were further selected from the 
21,026 unigenes by blast analysis with the e-value cut-off of 
1 x E~ s , to remove paralogous information and to build one- 
to-one relationships with the A. thaliana genes. It is noteworthy 
that there were 12,311 expressed genes common to the data 
sets from all three species, corresponding to 71% of the total 
expressed genes in A. thaliana, 64% in A. halleri and 87% in 
B. rapa, indicating that a large proportion of the expressed 
genes were common to papilla cells across the species boundary 
in Brassicaceae. Indeed, among the top 100 most highly 
expressed genes, 82, 82 and 78 genes in A. thaliana, A. halleri 
and B. rapa, respectively, were common to all three species 
(Supplementary Tables S1-S3). 

These gene sets were next sorted into six functional cate- 
gories, taking into consideration several aspects of the biolo- 
gical role, function and physiology of stigmas, i.e. signal 
transduction, cell wall-related, stress/defense, metabolism, ion 
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transport and transcription: this categorization is based on a 
functional classification of the pistil-expressed genes according 
to the GO annotations (Tung et al. 2005, Allen et al. 201 0). From 
the 17,240, 19,260 and 14,099 genes expressed in papilla cells 
of A. thaliana, A. halleri and B. rapa, 6,468, 7,481 and 5,700 
unigenes were classified into six groups according to their 
annotations and GO terms, respectively. General classification 
based on GO terms of the papilla-expressed gene data sets 
revealed that the proportion of categories showed broad simi- 
larities among species, in terms of functional classes (Fig. 6; 
Supplementary Table S6). These results reflect a degree of 
functional conservation of papilla cells between Arabidopsis 



A. thaliana 

18.8% 8.1% 



A. halleri 

^^Jki% 

1 9.9% 8.6% 



B. rapa 

^^05.9% 
1 9.4% 8.4% 



4.3% 



15.6% 



4.0% 



15.9% 



3.9% 



15.3% 



46.7% 



Signal transduction 
Metabolism 



45.5% 



Cell wall 
Ion transport 



47.0% 



Stress/defense 
Transcription 



Fig. 6 Classification of papilla-expressed genes by gene ontology, in 
A. thaliana, A. halleri and B. rapa. The pistil-expressed gene sets were 
sorted into six functional categories based on a functional classifica- 
tion of the genes, according to the gene ontology (GO) annotations 
(Tung et al. 2005, Allen et al. 2010): signal transduction, cell wall- 
related, stress/defense, metabolism, ion transport and transcription. 



and Brassica, and this is consistent with previous results in 
stigma (Tung et al. 2005, Allen et al. 2010). They also support 
the theory that mechanisms for stigmatic function are similar in 
different species, depending on the degree of species differen- 
tiation (Allen et al. 2010). A previous study indicated that the 
expression profiles of stigma-expressed genes in maize were 
more similar to those of rice than those of Arabidopsis 
(Swanson et al. 2004, Li et al. 2007, Xiao et al, 2012). 

In contrast to the above consideration of the similarity 
between the three species, 1,081, 2,376 and 581 papilla- 
expressed genes from A. thaliana, A. halleri and B. rapa, respect- 
ively, were detected in one of the species but not the other 
two (Fig. 5). This is presumably a result of species differenti- 
ation in Brassicaceae. 

Biological and functional characteristics of papilla 
cells in Brassicaceae 

To characterize the biological and functional roles of papilla 
cells in Brassicaceae, GO terms were assigned and sorted for 
the 12,311 expressed genes common to all three species. 
Among these, 4,309 distinct GO terms were annotated, con- 
sisting of 2,214 biological processes, 416 cellular components 
and 1,679 molecular functions (Supplementary Table S7). The 
10 most highly represented GO terms in each category of bio- 
logical process, cellular components and molecular functions 
are shown in Fig. 7. 

Metabolic process is the highest within biological processes 
and not the highest in all categories (Fig. 7). In 'metabolic 
process', highly represented GO terms included members 
of lipid metabolism. For example, among our data sets, 



Biological process (BP) 
Molecular function (MF) 
Cellular component (CC) 




Fig. 7 Functional categories of the 12,311 genes common to three Brassicaceae species. GO terms were assigned and sorted for 12,311 expressed 
genes common to three Brassicaceae species. 
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Extracellular Lipase 4 (EXL4) was expressed almost identically in 
papilla cells in all the species studied (Supplementary Table 
S7). EXL4 is a member of the GDSL lipases, which are active in 
hydrolysis and synthesis of lipids and/or esters (Updegraff et al. 
2009). GDSL lipases have multiple physiological roles such as 
development, flowering and seed germination in plants, and 
EXL4 is involved in the rapid initiation of pollen hydration. 
The ex/4 disrupted mutation results in slower pollen hydration 
on the stigma and decreased competitiveness in pollination, 
suggesting that changes in lipid composition at the pollen- 
papilla interface are essential for pollen hydration on papilla 
cells. Lipids have a broad range of biological functions, such as a 
component of the cell membrane and cuticle layer in epidermal 
cells, a mediator of signal transduction across extracellular and 
intracellular membranes, and an energy source for cellular func- 
tions (Wolters-Arts et al. 1998, Wolters-Arts et al. 2002), and 
the high proportion of lipid metabolism genes highlights their 
important role both in general cellular function (papilla cell 
development) and specific functions (pollen adhesion, pollen 
tube growth and pollen tube directional guidance and 
signaling). 

Another represented GO category within 'metabolic pro- 
cess' was associated with 'cell wall-related', including a 
member of the pectin methylesterase (PME), pectin lyase-like 
superfamily (PPME) and glycosyl hydrolase family (Bosch et al. 
2005) (Supplementary Table S7). The plant cell wall is com- 
posed of polysaccharides and proteins, and the polysaccharides, 
such as cellulose, hemicellulose and pectin, constitute the 



major part of the cell wall (Lehner et al. 2010, Komaki and 
Sugimoto. 2012). The varying proportions and interactions of 
the polysaccharides ensure the strong and flexible properties 
of cell walls (Daher et al. 2011). In Brassicaceae, the cell wall of 
papilla cells is composed of two layers, layer I (cuticle) and layer 
II (pectocellulosic wall) (Elleman et al. 1988). PPME and PME are 
proposed to mediate cell wall loosening and expansion, thus 
enabling pollen tube growth through the papilla cell wall 
(Bosch and Hepler 2005, Jiang et al. 2005, Renault et al. 2011). 
Application of an exogenous PME induces thickening of the 
apical cell wall and dissipation of the intracellular tip-focused 
Ca 2+ gradient in papilla cells (Hulskamp et al. 1995, Takamatsu 
and Takagi 2011). Thus, papilla cells in Brassicaceae have a 
broad range of metabolic turnover and cell wall remodeling 
systems, and genes in lipid metabolism and cell wall-related 
categories are expected to be important players in successful 
fertilization, through cell wall remodeling and cell-cell signaling 
interactions in papilla cells. 

The second category of biological process represented in the 
data set was 'regulation of transcription'. Based on the Database 
Arabidopsis Transcription Factors (DATF: http://datf.cbi.pku. 
edu.cn/), several families of transcription factors (TFs) were 
represented in the 12,311 common genes, and 754 TF genes 
were identified (Supplementary Table S8). Among these, 
six families of zinc finger (C2H2, C3H), PHD, MYB, bHLH and 
bZIP groups were represented (Fig. 8A). HECATE encodes a 
basic helix-loop-helix (bHLH) TF which controls the repro- 
ductive development of the transmitting tract and floral 



A 




Fig. 8 Distribution of TFs and receptor-like kinases (RLKs) in the 12,311 genes common to three Brassicaceae species. TFs were assigned by the 
Database Arabidopsis Transcription Factors (DATF: http://datf.cbi.pku.edu.cn/) for the 12,311 common genes. RLKs were identified and 
categorized among genes that have GO terms 'plasma membrane' within cellular components. 
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organ (Gremski et al. 2007). No Transmitting Tract (NTT) 
encodes a C2H2/C2HC zinc finger TF expressed in the transmit- 
ting tract, which is required for development of both the trans- 
mitting tract and the funiculus of ovules (Crawford et al. 2007). 
Other key genes were also identified: HALF FIELD (HAF) 
(Crawford et al. 2007, Crawford and Yanofsky 2011), 
Brassinosteroid enhanced expression 7 (BEE1) and BEE3 
(Friedrichsen et al. 2002). HAF, BEE1 and BEE3 encode closely 
related bHLH TFs that act redundantly to specify reproductive 
tract tissues, where they are expressed in distinct but overlap- 
ping patterns (Crawford and Yanofsky 2011). Thus, although 
little is known about the TF-controlled genetic network in re- 
productive development, the data sets of TF are expected to be 
a valuable resource for the study of genetic networks in papilla 
cell development in Brassicaceae. 

The predominant GO terms in the cellular components 
domain were 'chloroplast', 'plasma membrane' and 'cytosol'. 
A high proportion of the category 'chloroplast' is not likely to 
be papilla specific, because the proportion of the GO term 
'chloroplast' was uniformly high in other tissues such as leaf 
and stamen (data not shown). Several studies have shown that 
proteins integral to the plasma membrane, i.e. receptor-like 
kinases, play important roles in pollination, disease resistance, 
development and hormone response (Gaude and Dumas 1984, 
Hiscock and Allen 2008, Kato et al. 2010). The plant receptor- 
like kinases are transmembrane proteins that transduce exter- 
nal messages into the cell through their extracellular domains 
and intracellular kinase domain (Takasaki et al. 2000, Chen et al. 
2009, Shinya et al. 2010). In our papilla-expressed gene data sets, 
protein kinase families and leucine-rich repeat receptor-like 
protein kinases (LRR-RLKs) were the most highly represented 
in the 'plasma membrane', many of which also corresponded to 
'binding (protein binding, ATP binding, zinc ion binding, DNA 
binding)' in molecular function. Thus, it is likely that kinase 
families function in many aspects of papilla cells, such as co- 
ordination of pistil development, defense response and facili- 
tation of pollen-pistil communication, with downstream 
mechanisms in the cytosol. 

To understand further the role of protein kinase in papilla 
cell development and reproduction, we focus here on LRR-RLK 
in the papilla-expressed gene data set. Among 12,311 common 
genes, 220 LRR-RLKs were classified into 56 subfamilies (Fig. 8B; 
Supplementary Table S9). The first and fourth subfamilies 
included 22 cytoplasmic receptor-like kinases (RLCKs), which 
may modulate pollen tube growth and pollen-pistil inter- 
action. Botrytis-induced kinase 1 (BIK1) is a member of the 
RLCK Vila subfamily and BRI1 -associated receptor kinase med- 
iating brassinosteroid (BR) signaling (Clouse and Sasse 1998, 
Muschietti et al. 1998, Zuhua et al. 2000). The expression of 
BIK1 was highly represented in the common genes. 
Furthermore, the BR-responsive genes BEE1 and BEE3, down- 
stream factors of the BRI1 receptor complex, were also repre- 
sented. As mentioned above, these are TFs that regulate 
stigmatic cell development. Thus, BR signaling contributes to 
the pathway of papilla cell development, with a coordinated 



system of BIK1, BEE1, BEE3 and other as yet unidentified factors 
(Crawford and Yanofsky 2011). The second and third subfami- 
lies included genes encoding LRR-RLKs that may be involved 
in the regulation of pollen tube growth. In tomato, Solanum 
\y coper sicum, for example, pollen receptor kinases (PRKs) act in 
pollen-pistil interactions, mainly perceiving signals from pistil 
to regulate pollen tube growth (Muschietti et al. 1998, Tang 
et al. 2002, Johnson and Preuss 2003). Among the papilla cell- 
expressed genes, we identified several kinases, and these may 
play a role in coordinating pollination events and processes in 
pollen-papilla communication. Further analysis of this protein 
family will provide new insights into how pollen-papilla inter- 
actions are regulated and which are necessary for the signal 
cascade of pollination in Brassicaceae. 

Conclusion 

Pollination is an essential step for sexual reproduction, and 
its regulation is an important mechanism to maintain genetic 
barriers and diversity in plant species. Because correct function- 
ing of papilla cells is critical for successful reproduction in 
flowering plants, we focused on the establishment and charac- 
terization, on the basis of gene expression, of the molecular 
components of cellular events in papilla cells. From this 
approach, >60% of the expressed genes were found to be 
common to papilla cells of three Brassicaceae species, highlight- 
ing the functional commonalities of papilla cells of Brassicaceae 
in terms of development and reproduction. These mechanisms 
included, in particular, production of certain metabolites, 
TFs, receptor kinases and signaling cascades, which requires 
a precise gene network and molecular control to ensure cor- 
rect physiological, cellular, developmental and reproductive 
processes in papilla cells of Brassicaceae. Thus, papilla cells 
have complex cellular systems, which are regulated by a diver- 
sity of molecular players. In particular, LM-RNA-seq data reveal 
the wide variety of genes in each functional class that are 
involved in these processes. These data will be especially 
useful in studies of the spatial and temporal balance of gene 
expression that is critical for successful plant reproduction and 
development. 

The LM-RNA-seq data reported here highlight the effect- 
iveness of cell type-specific gene expression analysis, and these 
global gene expression data from papilla cells will be a valuable 
knowledge base for future research into plant reproduction, 
especially in pollination biology. This will shed light on the 
dissection of molecular mechanisms and components of 
pollination in papilla cells, aiding our understanding of the 
plant reproduction mechanism. 



Materials and Methods 




Plant materials 

Arabidopsis thaliana ecotype Columbia (Col-0) and its close 
relative A. halleri were grown in a growth chamber 
(BIOTRON LH-240S, NK system) at 22°C under 16 h light/8 h 
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dark photoperiod conditions. The individual of A. halleri subsp. 
gemmifera (Matsumura) O'Kane & Al-Shehbaz used in this 
study was collected from a population located in Tadaginzan, 
Inagawa, Osaka, Japan, by K.K.S. [population of halgem2 
described by Shimizu-lnatsugi et al. (2009) and Tsuchimatsu 
et al. (2012)]. It was subjected to five rounds of self-fertilization 
by bud pollination in the laboratory to reduce heterozygosity, 
and is the clone of the individual W302 with S-haplogroup A 
used by Tsuchimatsu et al. (2010). Brassica mpa accession 
Chiifu-401-42, a DH Chinese cabbage line, was grown in a green- 
house under normal conditions, at Chungnam National 
University, Daejeon, Korea. 

Tissue fixation and paraffin embedding 

Pre-chilled chemical fixatives (100% ethanol, 75% ethanol/25% 
acetate, 60% ethanol/40% acetate and 100% acetone) were 
used for fixation. Pistils were collected from floral buds before 
anthesis, and 3-4 mm tips including the stigma were picked 
and placed immediately in each fixative solution, with a 
vacuum on ice for 5 min twice. Tissue fixation and paraffin 
embedding were performed using a microwave processor 
(LabPulse H2850, Energy Beam Sciences), as described by 
Takahashi et al. (2010), with minor modifications. Each fixative 
was replaced with a fresh aliquot, and samples on ice were 
microwaved at 150 W, 37°C for 15 min three times. Fixed sam- 
ples were dehydrated by serial dilutions of 70, 80, 90, 100% and 
absolute ethanol, by microwaving at 250 W, 58°C for 90 s. In the 
case of samples in 100% ethanol and 100% acetone, the dehy- 
dration step was omitted and was replaced with 50% ethanol/ 
50% 2-propanol or 50% acetone/50% n-butanol, respectively, in 
infiltration steps. For paraffin embedding, each fixative was 
replaced by 50% paraffin/50% 2-propanol by microwaving at 
250 W 58°C for 10 min. In the case of samples fixed in 100% 
acetone, the fixative was replaced by 50% paraffin/50% n-buta- 
nol. All samples were then microwaved in 100% paraffin 
wax (Paraplast X-TRA, Fisher Scientific) at 150W, 58°C for 
10 min, followed by microwaving at 150W, 58° C for 30 min 
four times. The specimens were then placed on a plastic 
dish filled with melted paraffin wax. The embedded specimens 
were cooled to room temperature and the paraffin blocks were 
stored at 4°C 

Preparation of paraffin-embedded sections 
and LM 

The paraffin -em bedded samples were cut into 6-8 |im thick 
sections using a microtome (RV240, YamatoKoki). Serial paraf- 
fin-embedded sections were mounted on PEN membrane 
frame slides (Applied Biosystems, Life Technologies) with 
RNAsecure Reagent (Ambion, Life Technologies), dried on a 
heating plate at 57°C for 30 s to 1 min, and then further dried 
at 4°C for 1 h. To remove paraffin, slides were immersed in 
Histo-clear II (National Diagnostics) for 5 min at room tempera- 
ture twice, and dried completely at room temperature. LM was 
performed using an Arcturus XT Laser Capture Microdissection 
System (Applied Biosystems, Life Technologies). The procedure 



of LM is described in Suwabe et al. (2008). Briefly, selected 
areas were captured by an infrared laser onto CapSure Macro 
LCM Caps (Applied Biosystems, Life Technologies) and were 
subsequently cut by UV laser. The papilla cells that fused 
to the LCM cap were collected by removing the cap from the 
tissue section. 

RNA extraction, amplification and cDNA 
preparation 

Total RNAs of papilla cells collected by LM were extracted using 
the Pico-Pure RNA isolation kit (Applied Biosystems, Life 
Technologies). The quantity and quality of total RNA were as- 
sessed using an Agilent 2100 Bioanalyzer and RNA 6000 Pico kit 
(Agilent Technologies). The quality of the total RNA was judged 
by the RIN, which was calculated by 2100 Expert Software 
(Agilent Technologies). The RIN evaluates the quality of total 
RNA on a scale from 1 to 10, in which 10 corresponds to intact 
RNA and 1 corresponds to extremely degraded RNA. In plants, 
RNA with a quality value of more than RIN = 6 is acceptable for 
gene expression analysis (Takahashi et al. 2010). Thereafter, the 
mRNA was linear-amplified by a RiboAmp HS PLUS RNA amp- 
lification kit (Applied Biosystems, Life Technologies), according 
to the manufacturer's instruction, through 1.5 cycles. An aver- 
age of 5 ng of total RNA from each sample was used for linear 
amplification. cDNAs were synthesized from the linear-ampli- 
fied mRNA using the Superscript III First-Strand Synthesis 
system (Invitrogen, Life Technologies). The quantity and quality 
of cDNA were confirmed using Qubit fluorometer and Qubit 
dsDNA HS Assay kits (Invitrogen, Life Technologies). 

RNA-seq by SOLiD 5500x1 

Preparation of a fragment library for SOLiD sequencing was 
carried out following the manufacturer's instructions 
(Applied Biosystems, Life Technologies). Briefly, cDNAs were 
amplified by 10 cycles of PCR and a fragment library prepared 
using the 5500 SOLiD Fragment Library Core kit. Adaptor liga- 
tion was conducted using 5500 SOLiD Fragment Library 
Barcode Adaptors. The quality of the cDNA fragment library 
was confirmed with a fragment peak of 200-300 bp and the 
absence of unreacting adaptors, by an Agilent 2100 Bioanalyzer 
and Agilent High Sensitivity DNA kit. Approximately 
0.7-1.0 pmol of the cDNA fragment library of each papilla cell 
sample was applied to the SOLiD EZ Bead System, and then 
subjected to sequencing by the SOLiD 5500x1 system (Applied 
Biosystems, Life Technologies), with a 1/6 plate scale per 
sample. Because of the linear amplification of mRNA, SOLiD 
RNA-seq was not directional. 

Sequence data processing and gene ontology 
analysis for unigenes 

Processing of sequence data for assembly was conducted using 
LifeScope Genomic Analysis Software (Applied Biosystems, Life 
Technologies), with a default parameter. Sequence reads with 
a quality of >Q10 were assembled to reference genome 
sequences (TAIR10 for Arabidopsis; http://www.arabidopsis.org/, 
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BRAD for Brassica; http://brassicadb.org/brad/). The Arabidopsis 
unigenes were obtained from the TAIR and the B. rapa unigenes 
were searched using the BLASTX program against TAIR10, with a 
cut-off e-value of 1 x E~ s , and assigned GO terms. Examples of 
visualization of mapping reads were displayed by Integrative 
Genomics Viewer (IGV) 2.3, for three representative genes, 
AT1G09690, AT1G53130 and AT2G41430, which were used for 
in situ hybridization (Supplementary Figs. S2-S4). 

qRT-PCR 

Total RNAs were extracted from stigma and leaf samples ac- 
cording to the method described above. The concentration and 
purity of total RNAs were determined by spectrophotometric 
analysis. Total RNA was reverse-transcribed using the 
Superscript III First-Strand Synthesis system with oligo(dT) 
primer, according to the manufacturer's instructions. qRT- 
PCR was performed using the SYBR Green Real-Time PCR 
Master Mix (Toyobo) with a Bio-Rad CFX96 Real-Time 
Detection System. Quantitative validation in three replicates 
was calculated using the delta-delta threshold cycle relative 
quantification method. A gene coding for ubiquitin-conjugat- 
ing enzyme 21, UBC21, was used as a control. Sequence infor- 
mation of all primers is listed in Supplementary Table S10. 

In situ hybridization 

Pistils for in situ hybridization were collected 1 d before anthe- 
sis. Both antisense and sense probes were synthesized using a 
T7/SP6 digoxigenin (DIG) RNA labeling kit (Roche), according 
to the manufacturer's instructions. In situ hybridization was 
performed according to Masuko et al. (2006), and tissue sec- 
tions were observed under a light microscope by using SteREO 
Discovery V20 (Carl Zeiss, Jena, Germany). 



Supplementary data 

Supplementary data are available at PCP online. 
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